Classical dynamics of a two-species condensate driven by a quantum field 
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We present a stability analysis of an interacting two-species Bose-Einstein condensate driven by 
a quantized field in the semi-classical limit. Transitions from Rabi to Josephson dynamics are 
identified depending on both the inter-atomic interaction to field-condensate coupling ratio and the 
ratio between the total excitation number and the condensate size. The quantized field is found 
to produce asymmetric dynamics for symmetric initial conditions for both Rabi and Josephson 
oscillations. 
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I. INTRODUCTION 

Recently, a classical bifurcation at the transition from 
Rabi to Josephson dynamics has been observed experi- 
mentally in a rubidium spinor-Bosc-Einstcin condensate 
(BEC) driven by an electromagnetic field through a two- 
photon coupling process [l[ . A simplified model to access 
the collective dynamics of a similar full quantum system 
is that consisting of Nq interacting qubits driven by a 
quantum field through a one-photon process and given 
by the Hamiltonian, 



H = 5J^ 



XN-'^^a + a^)X, (1) 



where the detuning S is the difference between the BEC 
ground state hyperfine transition frequency w and the 
quantized field frequency ujf, S — uj — ujf. The orbital 
angular momentum representation, with the population 
difference and coherences given by Jz , Jx and Jy , in that 
order, has been chosen to describe the ensemble of dipole- 
dipole interacting qubits, while the field is described by 
the creation (annihilation) operator, a (d^). The field- 
ensemble coupling and the intra-ensemble coupling are 
given by the constants A and 77, respectively. 

A standing conjecture states that entanglement in non- 
linear bipartite systems can be associated with a fixed- 
point bifurcation in the classical dynamics [H, thus pro- 
viding a link between classical and quantum regimes. En- 
tanglement is a fundamental quantum phenomenon Q 
and resource Interactions are essential to gener- 

ate quantum correlations; 6.17., the Lipkin-Mcshkov-Glick 
(LMG) model, originally proposed in nuclear physics 
produces maximal pairwise entanglement of qubits at a 
quantum phase transition of its ground state and may 
describe the Josephson effect in a two- mode BEC 0] . 

The Hamiltonian in Eq.([T]), hereby called LMG-Dicke 
(LMGD), can be obtained from the Gross-Pitaevskii 
equation describing a two-species BEC interacting with 
a quantum field by following a derivation similar to 
that found in Rcf. Experimental realizations pro- 
viding an assorted range of tunable parameters for the 
LMGD model may include a two-hyperfine-structure- 
defined-modes BEC coupled to a quantum cavity field 
mode through a one microwave photon process; e.g. 



trapped hyperfine ground states of a Sodium BEC in- 
side a microwave cavity Q . Arrays of interacting super- 
conducting qubits coupled to the quantum field mode of 
a coplanar waveguide resonator may be considered with 
the limitation that ensemble sizes are small 

The ground state phase transition of the LMGD model 
has been studied in the thermodynamic limit, Nq — >■ 00, 
within the rotating wave approximation (RWA) [3] and 
in the quantum regime, using coherent states for both the 
field and ensemble, without the RWA [T^l. These results 
show the existence of a finite size first order quantum 
phase transition and a second order super-radiant phase 
transition. In the quantum regime, we have shown 
that maximal shared bipartite concurrence of the ensem- 
ble may be obtained for weak coupling, as the LMGD 
Hamiltonian in the limit A — )• becomes the LMG model. 

Here, we present a steady state analysis of the equa- 
tions of motion for the system in the large ensemble size 
limit to explore its collective dynamics. First, the rela- 
tion between the quantum and classical field drives are 
discussed comparing the Rabi and Josephson oscillations 
that appear in both cases. Then, we find that the sym- 
metry of Josephson dynamics is broken by the quantum 
field. Finally, an actual pitchfork bifurcation point is 
found in the regime where the intra-ensemble interac- 
tion is larger than the field-ensemble coupling, rj ^ X; 
maximal entanglement is found in this regime for the 
quantum analysis [ll| . Our results may provide a deeper 
understanding about the collective dynamics of interact- 
ing qubits and another example in favor of the aforemen- 
tioned conjecture relating classical and quantum regimes 
for nonlinear systems. 



II. MODEL 

Starting from the LMGD Hamiltonian defined in 
Eq.©, by considering a large ensemble, Nq ^ 1, it 
is possible to approximate the expectation values using 
the thermodynamic limit where the system is consid- 
ered to be in a separable state composed by a coher- 



ent field, 



and coherent atomic state, \0,ip)', i.e.. 
- y^e^^t il) « {Nq /2) cose, (J±) ft 



2 



{Nj2)smee^^f, and (iV) = iV « n - (A^,/2) cos 6*. By 
defining a fractional population difference, z = cos 6*, an 
excitation ratio parameter, k = 2N/Nq, and a total phase 
variable, $ = (j> + (p, the effective Hamiltonian in this 
mean-field approximation, under the RWA in units of 
hNqX/2, is given by the expression, 

H ^ (^A + ^^ z+ [2(1 - (fc _ z)] cos $, (2) 

where the dimcnsionless couplings ratio A = rj/X, detun- 
ning A — S/X, and the shorthand notation for the mean 
value of an operator O = (O) has been used. 

The mean field Hamiltonian in Eq. ([5]) describes a 
non-rigid non-linear pendulum and is equivalent to the 
case of a BEC in an asymmetric double well via a phase 
TT-shift and a restriction given by {k — z) = 1/2, cf. 
Eq.(5) in [l^ . Under these restrictions, plus symmetry 
of the double well, A = 0, the model describes a Bose- 
Josephson junction supporting coherent oscillations, i.e., 
Rabi dynamics for macroscopic self-trapping modes and 
Josephson dynamics below and above a critical tunnel- 
ing strength given by Ac = 1 [13]. Due to the phase 
shift with respect to the Bose- Josephson junction, the 
so-called plasma and tt oscillations U3|. alrcadv observed 
experimentally for bosonic BECs [l|, will exchange 
places in the studied system; i.e. plasma and tt oscilla- 
tions will be located at $ = tt and $ = 0, in that order. 

Notice that under resonant quantum driving, S = 
u! — Ljf = leading to A = 0, all deviations from the 
symmetric classical driven system are due to the quan- 
tized field; the term (fc — z)^/^ in Eq.® depends on the 
parameter k, i.e. the ratio between the conserved total 
number of field and atomic excitations in the system and 
the size of the condensate. 



III. FIXED POINTS 




Figure 1. (Color Online) Fixed points, with $ = 0, tt, 
defined by the parameters set given by the dimensionless de- 
tunning, couplings ratio, fractional population difference and 
excitation ratio, {A, A, z, k} in that order, for excitation ratio 
parameters fc = 0.1 (light blue) and k = 10 (dark red). 
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The fixed points of the model, Eq.®, arc found from 
the mean-field version of the quantum equations of mo- 
tion, up to a Nq/2 scale factor. 



^z = [2(1 - 2^) (fc _ z)] ^ ^.^ 
d , X l + 2fcz-3z2 

— $ = S + KZ 

at 



V^[(l-z2)(fc_2:)]l/2 



cos$. (3) 



The fixed points coincide with the critical points as i = 
dH/d^ and $ = dH/dz. Stationary states are found 
for the phase variable values $ = 0, tt and the excitation 
parameter value 



k = 



iz' 



1 (l-z2)|(A + Az)| 



2z 



4^2 



A + Az|± (A + Az)2-4z 



1/2 



(4) 



Notice that, in order to obtain a real excitation ratio, fc, 
the allowed fractional population difference is bounded 



Figure 2. Fractional population difference, z, as a function of 
the coupling ratio. A, for two excitation ratios, = 0.1 (grey) 
and fc = 10 (black) for fixed points Fq (solid) and (dotted) . 
Rabi oscillations occur around the fixed points marked from 
A to C, Josephson oscillations around fixed points D and F, 
finally, plasma oscillations around those marked G and H. 
The dashed branch for the latter excitation ratio, k = 10, 
shows the fixed points acting as separatrix of the two localized 
oscillations in the Josephson regime, in particular the fixed 
point E acts as a separatrix of oscillations around fixed points 
D and F. 



to the range z € [—1, zJ\ U 1], where z± = [2 — AA± 
2(1- AA)i/2]/A2 sets the condition A < 1/A. Requiring 
population inversion, Jz = Nq/2, that is z = 1, leads 
to fc = 1, i.e. N ^ Nq/2 which means there should be 
no excitations on the field, n = 0, without restriction on 
the couplings ratio A. In a similar way for = —Nq/2: 
z ^ -l,k = -I, N ^ -Nq/2, n^Q. Notice that z± = 
requires A or fc — > oo. 
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Figure [T] shows the fixed points, denoted as with 
$ = 0, TT, defined by the parameters set {A, A, z, k}, for 
excitation parameters k = 0.5 and k — 10, i.e. the total 
excitations number is smaller than the number of qubits 
and vice-versa, in that order. Figure [2] shows the plane 
defined by the resonant case, A = 0, from Fig. [TJ The 
different alphabetical markers in Fig. [5] pinpoint regions 
with different dynamics. The markers A and B are lo- 
cated in a region where there are only two fixed points 
per parameter set and these fulfill zp^ = zp^ . The next 
set of markers belongs to a region where two different 
sets of fixed points are identified; the one for a low ex- 
citation ratio where there is one fixed point for plasma 
and TT oscillations each, zp^ ^ zp^; the other for a large 
excitation ratio where there is only one fixed point for 
$ = TT, and three for $ = 0, i.e. a transition from Rabi 
to Josephson dynamics occurred for tt oscillations. 



IV. RESONANT QUANTUM DRIVING 
EXAMPLES, A = 0. 

Figure [3] shows an example where Rabi dynamics dom- 
inate, i.e. A < 1. The coupling ratio A = is chosen to 
explore the behaviour of the system when the inter-qubit 
coupling is negligible. The addition of the quantized field 
breaks the symmetry in the trajectories of the equivalent 
classical field driven case [l^. Rabi oscillations localize 
around the fixed point A in the Southern Hemisphere 
of the Bloch sphere defined by ($(i),z(i)); a sample of 
trajectories with initial conditions on $(i = 0) = 0, are 
shown in Fig. [Sljb) where two particular initial condi- 
tions, z = —0.5 and z ^ —0.9, are highlighted; a few of 
the trajectories show an unbounded phase. For a larger 
excitation ratio, fixed point B in Fig. Efa), Rabi oscilla- 
tions for symmetric initial conditions are still asymmetric 
as shown by the highlighted trajectories in dashed and 
solid lines portrayed in FigI3Jb) ; it is only in the limit of 
large excitation ratio, fc — >■ oo, that symmetric dynam- 
ics for symmetric initial conditions are recovered. Notice 
that unbounded phase modes disappear. Figures [2Jd- 
e) show the mean value of the mean-field Hamiltonian, 
Eq.(I3]), where it is possible to see that zp^ = zp^ for both 
cases of excitation ratios. 

On resonance, A = 0, for any given value of the cou- 
pling ratio. A, the allowed fractional population differ- 
ence is in the range z G [— 1, 0) U [4/A^, 1]; i.e. all of the 
stationary states are self-trapped states, z 7^ 0, unless 
A — > cxj. For the minimum self-trapped positive frac- 
tional population difference, = 4/A^, the excitation 
ratio is given by k{z^) = (16 -f A"')/8A^. Furthermore, 
setting dk/dz = delivers the cubic, with one real root, 



a) A(k,z)=0 



A 



^zl 



izi 



1=0, 



(5) 



that yields the set of critical parameters {fcc±(A, Zc), Zc] 
delimiting the Josephson dynamics region where there 
are localized oscillations around two fixed points with 
positive and negative fractional populations. Figure |3Ja) 
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Figure 3. (Color online) On resonance, A = 0, (a) Fractional 
population difference, z, as a function of the excitation ratio, 
fc, for a given coupling ratio, A = 0. Alphabetical mark- 
ers show fixed points for excitation ratios fc = 0.1 (A) and 
fc = 10 (B). (b-c) Trajectories for tt oscillations around the 
fixed points A and B for initial conditions including a survey 
of initial fractional population difference and a initial phase 
$(0) = 0; highlighted trajectories sharing the same tone cor- 
respond to symmetric initial conditions z(fi) = 0.9 (dotted 
gray), z(0) = 0.5 (solid black), z{Q) = -0.5 (dotted black), 
z — —0.9 (solid gray), (d-e) Mean value for the Hamiltonian 
showing the critical points defined by A and B. 



shows a graph of the fractional population difference, 
z, as a function of the excitation parameter, fc, for an 
on-resonance coupling ratio A = 6, where the criti- 
cal parameters {kc, Zc} are shown. For a low excita- 
tion parameter, fc = 0.1, coherent oscillations localize 
in the Southern hemisphere around the fixed point C, 
FiglHa); nrost of these Rabi oscillations are unbounded 
in phase, FigHlJb). Plasma oscillations localize around 
the fixed point G, which is the minimum of the mean 
energy, Eq.(I2]), Figl^d). For a larger excitation parame- 
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Figure 5. On resonance, A = 0, fractional population differ- 
ence, z (solid line), as a function of the excitation ratio, k, for 
a very large coupling ratio, A = 5000. A bifurcation point ap- 
pears at the critical excitation ratio kc+ ~ A^/2 = 1.25 x 10^. 
The separatrix for the Josephson regime is shown as a dashed 
line. 



of the self-trapped modes is asymmetric; symmetric dy- 
namics for symmetric initial conditions can be recovered 
only in the limit A — ^ oo. 

Figure ([5]) shows that for large coupling ratios, A ^ 1, 
i.e. 77 3> A, it is possible to locate a pitchfork bifurcation 
point, kc+ « This condition, 77 ^ A, relates to 

the phase space region vifhere maximal shared bipartite 
concurrence in the qubit ensemble may be obtained in 
the quantum treatment of this model . The difference 
comes from the large excitation parameter ratio arising 
in this semi-classical analysis, fcc+ 3> 1, i.e. Nq <C as 

« 3/A2 « 1. 



Figure 4. (Color online) On resonance, A = 0, (a) Frac- 
tional population difference, 2, as a function of the excita- 
tion ratio, fc, for a given coupling ratio, A = 6. Alphabet- 
ical markers show fixed points for excitation ratios A: = 0.1 
(C,G) and = 10 (D-F,H). (b-c) Trajectories for tt oscilla- 
tions around the fixed points C, for low excitation ratio, and 
D, for large excitation ratio, for initial conditions including 
a survey of initial fractional population difference and a ini- 
tial phase "1?(0) = 0; highlighted trajectories sharing the same 
color correspond to symmetric initial conditions z(0) = 0.9 
(dotted gray), z{Q) = 0.5 (solid black), z(0) = -0.5 (dot- 
ted black ), z = —0.9 (solid gray), (d-e) Mean value for the 
Hamiltonian showing the critical points defined by (C-H). 



ter, k — 10, an extra branch of fixed points appears and 
Josephson dynamics showing two self-trapped oscillation 
modes around fixed points D and F occur. Fig. IH^a). 
A sampler of the trajectories for initial conditions along 
<i>(t = 0) = 0, FiglDJc), shows that the fixed point E is 
the separatrix of the two sets of self-trapped oscillations. 
Figure |31[d) shows the mean energy for the system and 
how the fixed points correspond to local maxima, D and 
F, minima, H, and saddle, E, points. The localization 



V. DISCUSSION 

In the given examples, we have shown how the ratio 
between the total number of excitations of the system 
and the condensate size, k ~ 2N/Nq, limit the region of 
phase space accessible to both Rabi and Josephson oscil- 
lations. Only in the limit when the total number of field 
and atomic excitations in the system is infinitely larger 
than the size of the condensate, N ^ Nq/2 leading to 
k 00, the whole phase space is available for the oscil- 
lations. As the value of k diminishes, i.e. the number of 
field and atomic excitations become less, the available re- 
gion of phase space reduces, leading to more trajectories 
unbounded in phase for Rabi oscillations. 

The excitation ratio parameter, fc, may be highly tun- 
able as it relates to the population of the atomic species 
and the number of photons in the quantized driving field, 
all of them highly controllable in contemporaneous exper- 
imental realizations. 

More interesting is the asymmetric bifurcation of the 
nonlinear dynamics, with the displacement of the separa- 
trix from the typical value of z = 0, for values close and 
above the critical coupling ratio Ac = 1. A symmetric 
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pitchfork bifurcation is recovered for an infinitely large 
coupling ratio, A = 77/A — >■ 00; i.e., the LMGD Hamil- 
tonian, Eq.lJT]), becomes simply a LMG model, A 0. 
Hence the breaking of this particular nonlinear symme- 
try comes as a result of the quantized field driving. The 
behavior of the asymmetric bifurcation may be proved 
by varying the inter-atomic interaction, rj, with respect 
to the atomic-field coupling, A, with the condition that 
r] > X and the extra restrictions provided above to see 
Josephson dynamics. 

Notice that both the inter-atomic interaction and the 
atomic-field coupling are highly tunable parameters in 
current experimental realizations. Furthermore, in this 
particular scheme the weak coupling between a hyper- 
fine transition and a quantized microwave field within 
the RWA is desirable to explore large values of A. 

VI. CONCLUSION 

In summary, we have presented a stability analysis of 
the classical dynamics of a large ensemble of interacting 
qubits driven by a quantized field. There exists a transi- 
tion from Rabi to Josephson oscillations as in the classical 
field driven system. The quantized field produces local- 
ized asymmetric dynamics for symmetric initial condi- 
tions while the classical field produces symmetric dynam- 
ics for symmetric initial conditions. For low coupling ra- 



tios, depending on the excitation ratio, Rabi oscillations 
localize in the Southern Hemisphere of the correspond- 
ing Bloch sphere for low excitation ratios, and smoothly 
transit to cover both Hemispheres for large excitation 
ratios; completely identical coherent oscillations for sym- 
metric initial conditions are only recovered for infinitely 
large excitation ratios. For large coupling ratios, there 
exists a critical excitation ratio that defines two differ- 
ent sets of dynamics; below this critical excitation ratio, 
localized coherent oscillations, most of them showing an 
unbounded phase, appear in the Southern hemisphere of 
the Bloch sphere; above the critical excitation ratio, two 
localized, or self-trapped, oscillation modes characteristic 
of Josephson dynamics appear with asymmetric dynam- 
ics for symmetric initial conditions and a separatrix with 
a non-null fractional population difference. A pitchfork 
bifurcation, somehow similar to the classical field driven 
case, is found for large coupling and excitation ratios, in 
these cases the separatrix corresponds to ensemble states 
with negligible fractional population which are closer to 
their classical field driven equivalent. 
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